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It has been argued by Zurek and Kibble that the likelihood 
of producing defects in a continuous phase transition depends 
in a characteristic way on the quench rate. In this paper 
we discuss our experiment for measuring the Zurek-Kibble 
scaling exponent a for the production of fluxons in annular 
symmetric Josephson Tunnel Junctions. The predicted expo- 
nent is a = 0.25, and we find a — 0.27 ± 0.05. Further, there 
is agreement with the ZK prediction for the overall normali- 
sation. 

PACS Numbers : 11.27. +d. 05.70.Fh, ll.10.Wx, 67.40.Vs 



I. INTRODUCTION 

As the early universe cooled it is believed to have 
undergone a series of spontaneous phase transitions, 
whose inhomogeneities (monopoles, cosmic strings, do- 
main walls) have observable consequences, for structure 
formation in particular. These defects appear because 
the correlation length £ of the field (or fields) whose ex- 
pectation value is the order parameter is necessarily fi- 
nite for a transition that is implemented in a finite time, 
whether it be continuous or not. 

It is difficult to determine the microscopic dynamics of 
such fields but, using only simple causal arguments Kib- 
ble HH made estimates of this early field ordering, and 
the density of topological defects produced at Grand Uni- 
fied Theory (GUT) transitions at 10 _35 s. Unfortunately, 
because the nature of the field theories is not known with 
any reliability, and the effects of their evolution are not 
visible until the decoupling of the radiation and mat- 
ter 10 6 yrs later, it is impossible to provide unambiguous 
checks of these predictions. However, causality is such a 
fundamental notion that Zurek suggested that iden- 
tical causal arguments, with similar predictions, were ap- 
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plicable to condensed matter systems for which direct ex- 
periments on defects could be performed. In addition to 
their intrinsic interest for a better understanding of the 
dynamics of transitions in condensed matter, the hope 
is that successful tests of these predictions can lead to 
a better understanding of phase transitions in quantum 
fields. 

Several experiments in condensed matter systems have 
already been performed [p[~|Tl| to test the Zurek-Kibble 
predictions, with mixed results. It is these predictions 
that we have tested here, using annular Josephson tunnel 
junctions, for which the defects are fluxons. 

This paper is organised as follows. In the next section 
we discuss the Zurek-Kibble (ZK) scenario. Successive 
sections then give the predictions for annular Joseph- 
son Tunnel Junctions (AJTJs), show how the fluxons are 
measured, the nature of the samples, describe the experi- 
mental setup and, finally, present the measurements and 
their agreement with the predictions. An abbreviated 
description of the experiment has been given elsewhere 
12 1 . Details of the theoretical analysis can be found in 
13,^|, in which an earlier experiment [ jl5| by two of us 
(R.M and J.M) was analysed, to demonstrate its com- 
patibility with the ZK analysis, even though it had not 
been performed with this in mind. 



II. ZUREK-KIBBLE CAUSALITY 

Consider a system with critical temperature T c , cooled 
through that temperature so that, if T(t) is the temper- 
ature at time t, then T(0) = T c . T(0) = —T c /tq defines 
the quench time tq. 

There are several ways H of formulating the Zurek- 
Kibble causality bounds, but they all depend on the fact 
that, as the transition begins to be implemented, there is 
a maximum speed c(t) = c(T(t)) at which the system can 
become ordered. For relativistic quantum field theory 
(QFT) c is the speed of light. For superfluids, c(t) is the 
speed of second sound, vanishing at t = 0. For Josephson 
tunnel Junctions c(t), which depends on the nature of the 
junction, is the Swihart pa velocity. 
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Suppose that the equilibrium (adiabatic) correlation 
length £ad(t) = £ad(T(t)) diverges near t = as 



U(t) = £o 



t 



However, the true non-equilibrium correlation length £(t) 
can only change so much in a finite time, and does not 
diverge. Kibble and Zurek made two assumptions: 

• Firstly, the correlation length £ of the fields that 
characterizes the onset of order is the equilibrium 
correlation length £ = £ ffl d(t) a t some appropriate 
time t. 

• Secondly, we can measure £ experimentally by mea- 
suring the number of defects, assuming that the 
defect separation £^e/ = 0(£). 

There are several ways to estimate t (the 'causal time'), 
explicit or implicit in the early work of Zurek Q . Most 
simply, 

• £(t) cannot grow faster than c(t). This is true both 
before and after the transition. That is, t is defined 
by the condition that £ a d(t) ~ ~ c (t)- 

• The relaxation time for the long wavelength modes 
r(t) is defined by c(i) = £ a d(£)/ T W- From this 
viewpoint t is that time when we can return to an 
adiabatic regime, t w r(t). 

In simple systems these estimates agree, up to numeri- 
cal factors approximately unity Q. They give i of the 
form t ~ Tq 7 Tq , where ro <C tq is the cold relax- 
ation time of the longest wavelength modes, and the crit- 
ical exponent 7 depends upon the system. As a result, 
T Q 3> i S> To- Identifying the initial domain size and 
defect separation as in the second ZK assumption then 
gives 



£~£ad(*)=£o 



(1) 



where a = jv. This is very large on the scale of cold 
defects which shrink to size £ Q d(2/j n ) = O(£o), where 
Tfi n is the final temperature, and £0 is determined from 
the microscopic dynamics. We term a the Zurek-Kibble 
(ZK) scaling index. 

The analysis above is for large systems, of linear size 
L ^> £. For the experiment that we shall describe below, 
of Annular Josephson Tunnel Junctions (AJTJs) of cir- 
cumference C, we find that is C < £. In fact, we expect 
the conclusions to be equally valid for small systems, for 
which the relevant time might, incorrectly, seem to be the 
time when the coherence length becomes smaller than the 
system. The reason is that the causal bounds (0) are to 
be thought of as a short-hand for the underlying dynam- 
ics. At the microscopic level, causality along the lines 



above is not explicit, although encoded in the relevant 
dynamical equations. The picture is rather one of order 
being established through the growth of the amplitudes of 
long-wavelength instabilities. The earliest time at which 
we can identify defects from this viewpoint is when the 
order parameters have achieved their equilibrium mag- 
nitudes. Qualitatively, for simple models this time is in 
good agreement with the causal time t above. There is 
no real surprise in this. It has been shown by one of 
us [|l7],[l8| that, in general, the causal time and distance 
scales t and £ = £(£) are just as we would expect from 
dimensional analysis (in the mean-field approximation) 
and unstable modes grow exponentially, whereby the de- 
pendence of the causal time (and corresponding defect 
density) on the microscopic parameters is only logarith- 
mic. In the same way, the distance between field zeroes 
has the same scaling dependence on tq as £ of (Q) , up to 
logarithms. 

As for the production of defects, they are transitional 
regions between different system ground-states. For 
superfluids like 4 He and superconductors these transi- 
tional regions are flagged by zeroes of the scalar order- 
parameter field. For the case of JTJs we shall see 
that this is generalised to zeroes ( mod 1~k ) of the order- 
parameter field. Any field crossing zero (mod 2n) has the 
potential to mature into a defect. However, only when 
the transition is complete will the field configurations in 
the vicinity of the zeroes have the energy profile of a clas- 
sical defect, the solution to the classical field equations. 
Thus, before the causal time we now have a picture in 
which there is a fractal thermal fuzz of potential defects, 
whose density depends on the scale at which we look. 
By the causal time some of these have developed into the 
(scale-independent) defects that we see subsequently (see 
p7|,p8[ ). Because, from this viewpoint, we have many 
proto-defects jockeying to become the real thing, the rel- 
evant scale to compare to the system size is not £ but £0 , 
as before. Equivalently, counting field zeroes (mod 2-7T, 
or not) depends on the short-distance behaviour of the 
field correlations. We do not expect a problem as long as 
C ^> £0, as is always the case. 

We conclude by noting that this is a very different 
picture from that of domains freezing in as the transition 
is approached from above (t < 0), which is how causal 
bounds were originally posed ^,|) . What matters is that 
all these causal descriptions give results with the correct 
engineering dimensions. With this in mind, we keep the 
causal bound (^) as a convenient mnemonic. 



III. THE ZK PREDICTIONS FOR FLUXONS 

The order parameter of a Josephson tunnel junction 
at temperature T < T c is the phase difference (f> of the 
macroscopic superconducting quantum mechanical wave 
functions across the barrier. Using a Lagrangian formal- 
ism Gr0nbech- Jensen et al. [R3 have shown that, for an 
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annular JTJ with a distributed bias current 7, (j) obeys to 
the following perturbed Sine-Gordon equation (PSGE): 



d 2 4> 
dx 2 



1 d 2 (j) 
c^FjW 



1 



sm < 



: 7 + 



d(j> d 3 (/) 



'■(t) dt 







dx 2 dt 



(2) 



provided the width Ar of the annulus, of radius r, satis- 
fies Ar <C r and Ar <C A,/(T), the Josephson coherence 
length. In this case x measures the distance along the 
annulus, and c(T) is the Swihart velocity, a and /3 are 
the coefficients of the losses due to the tunnelling current 
and to the surface impedance, respectively. 

The boundary conditions for Eq.(||) are periodic pc| ] 
and derive from fluxoid quantization pi] ]: (f>(x + C) = 
4>(x) + 2irn, where C — 2nr is the circumference of the 
junction and the winding number n is an integer cor- 
responding to the algebraic sum of fluxons trapped in 
the junction barrier at the normal-superconducting (N- 
S) transition; n is a topological system constant, that 
is, only fluxon-antifluxon (FF) pairs can be created or 
annihilated as long as the junction remains in the super- 
conducting state. 

The classical fluxons are the 'kinks' of the Sine-Gordon 
theory. As with other models of defect formation, equa- 
tion (Q) is only valid once the transition is complete: 
therefore, we shall not use it to study the appearance of 
fluxons. However, it is sufficient to enable us, in the spirit 
of the Zurek-Kibble scenario, to identify Aj(T), diverging 
at T c , as the equilibrium correlation length £ a d{T) to be 
constrained by causality. Further, the Swihart velocity 
c(T) (with critical slowing down at T = T c ), measures 
the maximum speed at which the order parameter can 
change 

A detailed discussion of the ZK bounds has been given 
elsewhere by us |l3|,[l4]], and we refer the reader to these 
articles for more details. The JTJs in our experiment 
are symmetric, by which is meant that the electrodes 
are made of identical superconducting material with the 
same energy gaps and the same T c . For such JTJs [ p~3|Jl4| 
7 = v = 1/2. Therefore, at the time of their formation 

i = V T Q T 0: 

the separation of fluxons is expected, in the ZK picture, 
to be 



1/4 



(3) 



where £0 and To have the same meaning as before. It is 
the prediction (|j|) that will be tested in our experiment. 

In terms of the parameters of the JTJs, £0 has been 
inferred [Q as 



Co 



where J C (T) is the Josephson current density at temper- 
ature T. The parameter a is given in terms of the super- 
conductor gap energy and critical temperature and has a 
value between 3 and 5. If the thickness of the two super- 
conducting electrodes differs, the effective thickness d s is 
the harmonic mean of the individual thicknesses p2|| . As 
for to, it is given as to = £o/co, where Co defines the be- 
havior c(t) = co(t/TQ) 1 / 2 of the Swihart velocity for the 
system near T = T c . 



IV. MEASURING FLUXONS 

Once fluxons have appeared Eq.|| is relevant. A con- 
sequence of the periodic boundary conditions in AJTJs 
is that fluxons behave as relativistic particles on an in- 
finite lossy line. In the absence of any current through 
the barrier and/or externally applied magnetic field the 
fluxons experience a flat potential and therefore are in in- 
different equilibrium as far as the barrier is homogeneous 
and pin-hole free; in the reality, the barrier defects act 
as small pinning or repulsive potentials for the fluxons. 
Due to the losses, after the transient regime is over the 
fluxons are still. Unfortunately static fluxons are difficult 
to reveal, since according to the second Josephson law, 
any static phase profile 4>t(x,t) = does not alter the 
junction zero-voltage state that is also typical of a flat 
profile <p(x,t) — const., corresponding to the absence of 
any trapped fluxon. 

In contrast, whenever fluxons travel around an AJTJ 
they leave a clear signature on the junction current- 
voltage characteristic (IVC) and therefore are easily de- 
tectable. In fact, as soon as a bias current is fed to the 
AJTJ, the fluxons move as magnetic dipoles under the 
action of the resulting Lorentz force. The fluxon dynam- 
ics in long JTJs is a well known topic and has received a 
great deal of both theoretical and experimental attention 
in the last few decades. If the external bias is assumed to 
be uniform over the junction area, then, as a result of the 
balance between the externally supplied power and the 
internally dissipated power, the fluxons move with a con- 
stant speed: the larger is the external bias the greater is 
the fluxon speed, but never exceeding its relativistic limit 
sets by the Swihart velocity. The motion direction de- 
pends both on the current sign and fluxon polarity (i.e., 
whether n = ±1) for a given bias current. Fluxons hav- 
ing different polarities travel in opposite directions and 
are likely to annihilate when they collide at low speed. 

Quantitatively, if a fluxon travels around an AJTJ hav- 
ing a mean circumference C with a constant speed i>, 
then it has a angular speed uj = 2ttv/C and the phase 
(j> advances of 2ir each period T — 2it/uj = C/v. There- 
fore, according to the second Josephson equation, an av- 
erage voltage V develops across the junction equal to 
V = (d(f>/dt) = <S> /T = $ w/C that is proportional 
to the fluxon velocity. In other words, the presence of a 
travelling fluxon sets the junction in a finite voltage state 
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than can be easily measured on its IVC. 

By also changing the bias current through the barrier 
the voltage drop changes and a new branch called Zero 
Field Step (ZFS) appears on the junction IVC; the ZFS 
represents the relationship between the applied Lorentz 
force (proportional to the bias current) and the fluxon 
speed (proportional to the voltage). When N fiuxons 
travel around an AJTJ, the last expression is easily gen- 
eralized to give a junction voltage V — NQqC/v. In the 
last expression TV is the total number of travelling fiuxons 
and can be larger than the winding number n if FF pairs 
are travelling around the annulus. Therefore, we count 
the number of travelling fiuxons by simply measuring the 
voltage across the AJTJ. 

These properties makes AJTJs very competitive with 
respect to other solid state systems proposed to test the 
Zurek-Kibble mechanism. Our idea is to perform a large 
number of N-S transitions on the same AJTJ with no 
external current or magnetic field; at the end of each cycle 
the possible spontaneously generated fiuxons are static. 
Then we supply an external current that sets the fiuxons 
(if any) in motion around the annulus and measure the 
number of travelling fiuxons by a careful inspection of 
the junction IVC. Due to the annihilation of a fluxon- 
antifluxon pair, this idea works well as long as the chances 
to spontaneously generate two fiuxons are small. 

Fig. la, Fig. lb and Fig.lc represent the IVC of the same 
AJTJ with no fluxon trapped, with one fluxon trapped 
and with two fiuxons trapped, respectively. We note that 
with no trapped fiuxons the zero voltage current is very 
large and only FF current steps appear at finite voltage. 
In the other two cases the supercurrent is rather small 
(theoretically it should be vanishingly small in ideal, pin- 
hole free barriers) and large current branches can be ob- 
served at finite voltages corresponding to the fiuxons and, 
possibly, FF pairs travelling around the junction. 

V. THE SAMPLES 

High quality Nb/Al — Al ox /Nb JTJs were fabricated on 

0. 5 mm thick silicon substrates using the trilayer tech- 
nique in which the junction is realized in the window 
opened in a SiO insulator layer. Details of the fab- 
rication process can be found in Rcf. j23|. On each 
15 x 24 mm 2 chip four JTJs were integrated, of which 
three ring shaped junctions having a mean circumference 
C = 500 [im and a width Ar = 4 fim and one 4 x 500 /im 2 
overlap- type linear junction. The so called "idle region", 

1. e. the overlapping of the wiring layer onto the base elec- 
trode was about 3 pan for all the junctions. The thick- 
nesses of the base, top and wiring layer were 200, 80 and 
400 nm, respectively. 

For all samples the high quality has been inferred by 
a measure of the I-V characteristic at T — 4.2 K . In 
fact, the subgap current I sg at 2mV was small com- 
pared to the current rise AI g in the quasiparticle cur- 



rent at the gap voltage V g , typically AI g > 35I sg , the 
gap voltage was as large as V g = 2.76 ml/ and the maxi- 
mum critical current I c was larger than 0.55A/ g for the 
overlap type junction. Furthermore, the application of 
a strong enough external magnetic field in the barrier 
plane completely suppressed any Josephson structures in- 
dicating the absence of electrical shorts in the barrier. It 
is important to mention that i) no logarithmic singular- 
ity has been observed in the IVCs at low voltages and 
ii) the temperature dependence of the critical current 
was linear as the temperature T approached the criti- 
cal temperature Tq] both these observations assure us 
that the junctions are symmetric, i.e. no detectable dif- 
ference can be assumed between both the energy gaps 
Ai,2 and the critical temperatures Tq\,i of the junction 
electrodes 1 and 2 in the proximity of the barrier. The 
maximum Josephson current density Jj was of the order 
of lk A/ cm 2 corresponding to a specific barrier normal 
resistance of about 200 il/im 2 . 

Many samples have been measured. For clarity only 
two will be discussed here. The geometrical and elec- 
trical (at 4.2 K) parameters of the two selected annular 
junctions on different wafers are listed in Table I. They 
have the same geometry (both the base and top electrode 
have a hole concentric to the ring,), but differ in the crit- 
ical current densities, i.e., in the normalized mean cir- 
cumferences C/Xj. The critical current density has been 
calculated from the measured quasiparticle current step, 
AI g , at the gap voltage. The values of the barrier mag- 
netic thickness A = 180 nm has been used for numeric 
calculations. On each same chip a linear overlap- type 
junction with the same width, length and idle region was 
used in order to measure the junction Swihart velocity Cq 
with a geometry in which the effects of the self field are 
minimized. The value of cq — 1.4 x 10 7 m/sec, due to the 
effect of the idle region, is 1.5 times larger than that ex- 
pected for a bare junction. This value of cq corresponds 
to a value of 0.08F/m 2 for the barrier effective specific 
capacitance. 

The data in Table I show that both samples are high- 
quality, long (C ^> A,/) annular JTJs. 



Sample 


A 


B 


Mean circumference C'(fim) 


500 


500 


Width Ar{p,m) 


4 


4 


Zero field critical current I a (mA) 


33 


2.5 


Maximum critical current / m ax(m^4) 


39 


2.7 


Gap quasiparticle current step AI g {mA) 


88 


5.2 


-^max/ AIg 


0.45 


0.52 


Critical current density J C (A/ cm z ) 


3050 


180 


Josephson length Xj(fim) 


6.9 


28 


Normalized mean circumference C/Xj 


72 


18 


Quality factor V m (mV) 


49 


63 


Normal resistance i?jv(m£l) 


36 


610 


ZFS1 asymptotic voltage (pV) 


51 


53 



Table I 



4 



VI. THE EXPERIMENTAL SETUP 

In order to vary the quenching time in the broadest 
possible range, we have realized the experimental setup 
schematically shown in Fig. 2. 

A massive Cu block held to the sample holder by two 
thin Cu arms was used to increase the system thermal 
capacity. The chip was mounted on one side of this block 
and thermally loosely coupled to it by means of a 1 mm 
thick teflon sheet. On the other side of the Cu block, a 
thermoblock consisting of a 50 O carbon resistor and two 
thermometers in order to measure and to, if necessary, 
stabilize the Cu block temperature, was mounted in good 
thermal contact. Finally a small sized 100 O resistor, 
more precisely a surface mount resistor (SMR), was kept 
in good thermal contact with the chip by means of a 
small amount of vacuum grease. 

This system, due to the two heating elements placed in 
tight and loose thermal contact with the chip and with 
the proper choice of the thermal loads, allowed us to 
perform the sample quenching over two quite different 
time scales. In fact, by means of the resistor in the ther- 
moblock, a long time scale was achieved by heating the 
chip through the Cu block and the teflon sheet; on the 
contrary, a short current pulse through the surface mount 
resistor on the chip, attained much short thermal cycles. 

These two completely different quenching techniques 
provide timescale ranges that do not overlap, leaving 
a gap between 0.2 and Is, that would require a third 
quenching technique to be filled. We stress that just 
using a single sample holder with smaller heat capac- 
ity would not give us access to the same time-scales. We 
could not get quench times as short as those from the 
SMR and, even with a very small exchange gas pressure, 
could not get times as long as those with the larger sam- 
ple holder on using the mechanical pumps available to 
us. 

The whole system was kept in a vacuum tight can im- 
mersed in the LHe bath at He gas pipeline pressure. 
The pressure of He gas inside the can could be varied in 
order to modify the heat exchange between the chip and 
the environment and, in turn, the speed of the sample 
cooling. A solenoid was wound around the can to pro- 
vide a strong vertical magnetic field and Hclmholtz coils 
were instead placed inside the can to generate a weak 
horizontal magnetic field in the barrier plane in order to 
tune the critical currents of the annular junctions to their 
maximum values. 

The temperature dependence of the junction gap volt- 
age was exploited to monitor the temperature of the junc- 
tion itself during the thermal cycle. Fig. 3a and Fig. 3b 
show the digitally measured V g (t) for sample A current 
biased on the quasiparticle curve at 17.7mA, (i.e., at 
about one fifth of AI g at T = 4.2 K), for a slow and 
fast thermal cycle, respectively. In the case of Fig. 3a, 
a 100 mA current was fed to the thermoblock heater for 
about 4s in order to increase the junction temperature 



up to its critical value where the sample IVC becomes a 
straight line with a slope corresponding to the junction 
normal resistance Rn- In the case of Fig. 3b a 20 V high 
and 4tos wide, voltage pulse was applied at the SMR. In 
the Figs. 3 the time origin is arbitrarily set at the instant 
we began to feed the heating elements. It is important 
to observe that the time scale in Fig. 3a is about 50 times 
larger than that of Fig. 3b, although the curve shapes are 
quite similar. 

In our samples the current jump at the gap voltage is 
very steep and, at T — 4.2 K, the voltage changes by 
less than 1 % when the current is changed from 10 to 
90% of the total current jump Ai g and by less of 10% 
at T = 8.5 if. Furthermore, at T = 4.2 and 8.5 K and 
for this bias current, the junction voltage was equal to 
2.74 and 1.0 mV, respectively. Therefore, assuming that 
the electrode gap energies are equal Ai(T) = A 2 (T) = 
A(T), and that in the 4.2—8.5 K range we can neglect the 
thermal gap smearing, the analytical expression found 
by Thouless ]24| ] for the gap energy in a strong-coupling 
superconductor 



A(T) 
A(0) 



tanh 



A(T) T c 
A(0) T ' 



(4) 



also applies to the junction gap voltage that is propor- 
tional to it. The dashed lines in Figs. 3 indicate the volt- 
age threshold above which Eq^ can be used to relate the 
junction voltage at its temperature, without any signifi- 
cant loss of accuracy. We like to mention that an experi- 
mental proof of Eq.[| in Nb/Nb tunnel junctions was first 
evidenced by Broom 

Fig. 4a and Fig. 4b show the data reported in Fig. 3a and 
Fig. 3b, respectively, transformed according to Eq.||, as- 
suming for A(0) and T c the values 2.85 meV and 8.95 K, 
respectively, as found by Monaco et ol. p6| on similar 
JTJs. Now the dashed horizontal lines indicate the tem- 
perature threshold below which the temperature time de- 
pendence can be reliably accounted for by our measured 
data. The large noise at low temperature is the result 
of an amplification effect of Eq.||, according to which the 
temperature variation corresponding to a given energy 
gap variation becomes larger and larger as the tempera- 
ture becomes smaller. It would be very complicated to 
write the proper boundary conditions for the heat diffu- 
sion equation that would correctly model the full time 
dependence of the junction temperature. However, for 
our purposes we are only interested to the cooling pro- 
cess, and we successfully fit our data by a simple thermal 
relaxation equation: 



T(t) = Tf in + (T ln - T fm ) exp 



t-t 



(5) 



with only two fitting parameters to and r, being Tj„ and 
Tfi n fixed at 8.95 and 4.15 if, respectively. In Eq.[| to is 
the time at which T = Ti„ = T c and r is the relaxation 
time which sets the cooling time scale. The fitting curves 
are shown by the thick dashed lines in Fig. 4a and Fig. 4b, 
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and correspond to a thermal relaxation time r equal to 
3.6 and 0.073 s, respectively. The quenching time tq can 
be obtained from its definition: 



Tc 



dT 



\T=T C 



(6) 



giving tq = tTc / (Ti n — Tfi n ). For a thermal relaxation 
from the junction critical temperature down to the he- 
lium bath temperature we get tq ~ 1.9 r. Eventually, 
Fig. 5 display the values of the quenching times obtained 
with the process described above, both for the fast (black 
squares with right vertical scale) and slow (red circles 
with left vertical scale) cooling processes and for differ- 
ent values of the He pressure inside the can. We observe 
that, by changing the exchange gas pressure and using 
the two techniques, the quenching time can be changed 
over a quite large range starting from tenths to tens of 
seconds. At the end of this paper we will discuss how it is 
possible to extend this range in both directions; however, 
as we will show in next section, this range has shown to 
be large enough for our purposes. 



VII. THE MEASUREMENTS 

Quenching experiments were carried out in a double fi- 
metal shielded cryostat and the transitions from the nor- 
mal to the superconducting states were performed with 
no current flowing in the heaters and the thermometers. 
Both the junction voltage and current leads were shorted 
during all the thermal cycle. Furthermore, the heat sup- 
plied to the sample was such that the maximum temper- 
ature reached by the junction was made slightly larger 
than its critical temperature, say at about 10 K, in order 
to make sure that also the bulk electrode critical tem- 
perature (Tc — was overcome. In this case, ac- 
cording to Eqs.|| and n, the value of the quenching time 
results in a correspondingly smaller tq — 1.7 t. Due to 
the approximation made by using Eq.g and to the exper- 
imental uncertainty in the knowledge of the maximum 
temperature during each thermal cycle, the value of the 
quenching times has been determined with an overall ac- 
curacy as large as 5%. For each value of the quenching 
time, in order to estimate the trapping probability, we 
have carried out a set of 300 thermal cycles and at the 
end of each cycle the junction IVC was inspected in order 
to ascertain the possible spontaneous trapping of one or 
more fluxons. 

As we shall see later, the AJTJs are such that the 
ZK causal length £ > C by an order of magnitude when 
tq = Is. Increasing and decreasing tq by an order of 
magnitude changes £ by less than a factor of two. Thus 
the probability of finding a single fluxon after a quench 
is small. In the following we will focus our attention only 
on the probability P\ to trap just one fluxon, although 
a few times we found clear evidence of two and, more 
seldom, three homopolar fluxons spontaneously trapped 



during the N-S transition. However, these events were 
too rare to be statistically significant. 

Experimentally, we define P\ as the ratio between the 
number of times in which at the end of the thermal cy- 
cle the junction IVC looks like that shown in Fig. lb, i.e. 
with a tiny critical current and a large first ZFS, and the 
number of attempts. It is worth to mention here that 
in the case of simultaneous trapping of a fluxon and an 
antifluxon, they would annihilate and leave no track of 
their formation. Therefore, our definition of P\ is not rig- 
orous, but it is reasonable as far as the chances to trap 
two (homo or heteropolar) fluxons are negligibly small. 
Similarly, an IVC similar to that shown in Fig. lb could 
be the result of the simultaneous trapping of two flux- 
ons and one antifluxon, or the other way around, but 
this event is less likely in our experimental situation. For 
the sake of completeness, it must be added that, in some 
cases, the IVC displayed either a depressed critical cur- 
rent without ZFS structures or a ZFS with an enhanced 
critical current. We explain them as due to the trap- 
ping of Abrikosov vortices in the junction electrodes and 
nearby the barrier and we did not take in any account 
the occurrence of such events, since it is not known if 
and how the vicinity of Abrikosov vortices influences the 
fluxon formation. 



VIII. THE RESULTS 

When £ > C we estimate the probability of finding a 
fluxon in a single quench to be 



Pi 



i Co V ro 



(7) 



where, from Eq.(|J), a = 0.25. 

Fig. 6 shows on a log-log plot the measured probability 
P\ of a single fluxon trapping obtained by quenching the 
sample A 300 times for each value of the quenching time 
tq changed by varying the exchange gas pressure and 
by using both the fast and slow quenching techniques. 
We observe that the points are quite scattered mean- 
ing that the data are statistically poor. Further, for the 
reasons given earlier there is a gap between fast and slow 
quenches. Nonetheless, we have clear evidence that i) the 
trapping of a fluxon occurs on a purely statistical basis 
being identical the conditions for each thermal cycle and 
ii) the probability to trap one fluxon is larger when the 
transition is performed at a faster speed (smaller quench- 
ing time) in a qualitative accordance with the causality 
principle. More precisely, we can distinguish the point 
to the lower right of the graph corresponding to the 
slow cycle process which gives in the average a proba- 
bility to spontaneously trap a fluxon once every fourteen 
attempts, and the cloud of data to the upper left cor- 
responding to the impulsive junction heating and giv- 
ing an average probability of one successful event every 



G 



about six attempts. This suggests that possible temper- 
ature gradients induced by the SMR are not an impor- 
tant source of systematic error, since such gradients, with 
their slowly moving profiles, have a tendency to reduce 
defect production |27[ . 

Regardless of the data spread, as suggested if Eq.0 
holds true, we attempted to fit the data with an allomet- 
ric function P\ = a ■ t^ with a and b being free fitting 
parameters. We found that the best fitting curve, shown 
by the solid line in Fig. 6, has a slope b = 0.27±0.05. Such 
a value of b, although affected by a 20% uncertainty, is in 
good agreement with the forth root square dependence 
expected for symmetric junction. 

For the coefficient a we found the best fitting value 
of 0.1 ± 10% (tq in seconds). This is to be compared 

1 /4 

with the predicted value of Ct ' /£ • Sample A had a 
circumference C — 500 fim. Its effective superconductor 
thickness was d s « 250 nm. At the final temperature 
Tfin — 4.2 if, the critical current density was J c {Tfin) — 
3050 A/cm 2 and the Josephson length was Xj(Tfi n ) — 
6.9 fim. From this, and cq given earlier, we infer that £o ~ 
3.8 (Ltm and To ~ 0A7ps. This then gives Ct^ 4 /£o 
0.08 s 1 / 4 , in good agreement with the experimental value 
of b, given the fact that we only expect agreement in 
overall normalization to somewhat better than an order 
of magnitude level. At this level, such a result is immune 
to systematic error in one or other of the measurement 
processes. After the problems (discussed below) of the 
experiments discussed in to find (reliable) defects at 
expected densities, if at all, our experiment shows that 
the ZK estimate remains sensible. 

Similar measurements have been carried out for sam- 
ple B. Although not in contradiction with (jl|), the results 
were affected by a data scattering even larger than that 
found for sample A (shown in Fig. 6). This is due to a 
much smaller normalized length which, according to ([?]), 
translates in a expected probability P\ , for a given tq, 

about four times smaller (since Ct 1/4 /?o m 0.02 s 1 ' 4 for 
this sample) , far too small too get statistically significant 
data in reasonable times considering that these measure- 
ments are both very time and LHe consuming. In order 
to have data comparable with those of sample A, sample 
B would have required a quenching times 4 4 times larger. 
However the roughly measured probability P\ of 1 fluxon 
every 50-100 attempts is in fairly good agreement with 
the expected value. This shows that our fluxons are not 
spurious byproducts of the measurement mechanism. 

IX. COMMENTS, FUTURE EXPERIMENTS AND 
CONCLUSIONS 

We consider this experiment to give a strong confir- 
mation of the Zurek-Kibble predictions. We said in the 
introductory section of this paper that condensed matter 
experiments to test the ZK predictions had given mixed 



results, and it is interesting to put this experiment in 
that context. 

Prior to our experiment, seven other experiments had 
been performed to test (Q), five with fixed tq two 
with variable tq plO]. [ In addition, the experiment 
cited earlier Jl5| on JTJs by two of us (R.M and J.M) was 
compatible with (J3|) , although it had not been performed 
with a test of ([|) in mind. It was this that motivated the 
experiment described here.] 

Of those experiments with fixed tq, two were experi- 
ments (5]J^] on supcrfluid 3 He — B, which rely on the fact 
that, when it is bombarded with slow neutrons, energy 
is released which leads to a hot spot, with temperature 
T > T c , in the superfluid which then cools below T c . This 
leaves behind a tangle of vortices, the topological defects 
in this system, whose density can be measured. Since 
tq is fixed by the nuclear process, it is not possible to 
confirm the predicted value a — 1/4. However, with only 
a single data point conflating both normalisation and a 
both experiments are highly compatible with (Q). 

The remaining experiments with fixed tq were two 
|?],H on superfluid 4 He, and one || on high temperature 
superconductors (HTSC). 

In principle, the 4 He experiments MM, which use a 
pressure quench with a varying timescale tq to imple- 
ment the transition, could have allowed for a more com- 
plete test, in this case to confirm a — 1/3 (after renormal- 
isation group rescaling). Yet again, vortices are the rele- 
vant defects. In practice, the most reliable published ex- 
periment H sees no vortices. In this context, the vortices 
seen in an earlier 4 He experiment , at levels compatible 
with (|l|) , were most likely an artefact of the experimental 
setup. Further experiments on 4 He are underway. 

The fifth experiment §, on HTSC, measures total flux 
through a surface i.e. the variance in the topological 
charge, carried in this case by the Abrikosov vortices. 
The vortex separation of (Q) can be converted into a pre- 
diction for the variance, but no flux is seen in contra- 
diction with this prediction, despite the phase separation 
that is a prerequisite for the result being seen elsewhere 
p8|. There is no obvious explanation of this null result. 
An attempt to take gauge fields into account fully |2|] 
shows that there is an additional mechanism for vortex 
production in the thermal fluctuations of the magnetic 
field but, as yet, this seems insufficient to explain the 
result. Such a mechanism will not apply to the JTJs 
considered in our experiment. 

These early experiments have either provided one data 
point for (0 ), or have been null for whatever rea- 
son. Two subsequent experiments have permitted vary- 
ing quench rates and so an estimate for a. The most 
recent [[Tof involves the Benard-Marangoni conduction- 
convection transition, in which a homogeneous conduc- 
tion state is broken into an hexagonal array of convec- 
tion lines on heating. The defects here are not associated 
with the line zeroes of an order parameter field, and the 
viscosity-dependent a does not match the ZK prediction, 
most likely for that reason. The more relevant exper- 



7 



iment JO]] is carried out in a non-linear optical system, 
with complex beam-phase the order parameter, satisfying 
a time-dependent Ginzburg-Landau equation with drift. 
There has been much numerical analysis |]30f of time- 
dependent Ginzburg-Landau systems, which show agree- 
ment with the ZK predictions for scaling exponents. The 
control parameter in [ pT[ is not the temperature, but 
the light intensity. Increasing it leads to pattern for- 
mation (defects) at a critical value. The predicted scal- 
ing parameter a = 1/4 is recovered to good accuracy as 
&exp = 0.25 ± 0.02, but agreement with normalisation is 
not stated. 

Given this relatively poor success rate in confirming 
(|l|) we are considering a further experiment to measure 
the ZK scaling exponent, this time with manifestly non- 
symmetric AJTJs. In Jl3| and (l4| we observed that it 
should have been difficult to make JTJs truly symmet- 
ric, as those used here. However, in |l3) and 14 we had 
not appreciated how the difference between symmetric 
and marginally non-symmetric JTJs is smeared by the 
proximity effect of the Al within the insulating layer. 
Significantly non-symmetric JTJs require different fabri- 
cation techniques, but the value of a inferred from the 
same causal arguments is a = 1/1 ', very different from 
the value of 1/4 that we tested above. The data from 
our experiment is incompatible with a = 1/7. This does 
suggest that a further experiment, with markedly non- 
symmetric JTJs, should be performed. 

Our experiments have demonstrated that quenching 
time of the order of 1 second give a rather large proba- 
bility to trap one fluxons on AJTJs having a very large 
normalized length. However, very long junctions mean 
very large critical current densities that, in turn, require 
Josephson barriers so thin that their quality and unifor- 
mity is often spoiled; furthermore, in most cases, applica- 
tions require intermediate length junction or even small 
junction. For these reasons, it would be highly desir- 
able to compensate the reduced junction length with an 
increased quenching rate, as it is suggested by the find- 
ings for sample B. Therefore, we like to conclude the 
paper with some comments on the possible technical im- 
provements that would allow to test the ZK predictions 
over a broader quenching time range. Firstly, tq can 
be trivially heightened by increasing the Cu block ther- 
mal capacitance. On the contrary in order to lower the 
quenching time, that is to make the N — S transition 
faster, it is needed to resort to new techniques since the 
maximum power that can be dissipated by the surface 
mount resistors sets an obvious lower threshold on tq. 
One possible way to reach this goal is to perform the 

junction thermal cycle by means of light pulses. Light 
dissipates inside the superconducting electrodes, but not 
in the substrate providing a local junction heating that 
will relax much faster to the background temperature. 
We estimate that, by using a properly focussed pulsed 
light beam, the quenching time scale can be reduced to 
the ixs range. 
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X. TABLE AND FIGURE CAPTIONS 

Table I Geometrical and electrical parameters of two se- 
lected annular Josephson tunnel junctions at 4.2 
K. 

Fig. 1 Low voltage part of the experimentally measured 
current-voltage characteristics of the same annu- 
lar Josephson tunnel junction a) without trapped 
fluxons, b) with one trapped fluxon, and c) with 
two trapped fluxons. For each current branch the 
corresponding number of travelling fluxons F and 
fluxon-antifluxon pairs FF is indicated. 

Fig. 2 Sketch (dimensions are not to scale) of the cryo- 
genic insert developed to perform the junction ther- 
mal cycles with a time scale changing over a broad 
range. 

Fig. 3 Digitally measured time dependence of the junction 
gap voltage during the a) "slow" and b) "fast" ther- 
mal cycle. For these measurements the junction 
was biased at about 1/5 of the gap quasiparticle 
current step. The horizontal dashed lines indicate 
the voltage threshold above which Eq.|| can be used 
to relate the junction voltage at its temperature. 

Fig. 4 Time dependence of the junction temperature dur- 
ing the a) "slow" and b) "fast" thermal cycle ob- 
tained from the data of Figs. 3 transformed accord- 
ing to Eq.^|. The horizontal dashed lines indicate 
the temperature threshold below which the temper- 
ature time dependence can be reliably accounted 
for by our measured data. Furthermore, the thick 
dashed lines are best fitting curve of the cooling 
process according to the thermal relaxation expres- 
sion EqJ|. 

Fig. 5 Quenching time tq as a function of the He gas 
pressure inside the vacuum tight can. By changing 
the pressure of the exchange gas inside the can, 
the system thermal constants are varied. The solid 
squares refer to the left vertical scale while the solid 
circles refers to the right vertical scale. 

Fig. 6 Log- log plot of the measured probability Pi to trap 
one fluxon versus the quenching time tq . The solid 



line is the best fitting curve found assuming a power 
law dependence as suggested by Eq.JjJ To a good 
degree of approximation the fit is in agreement with 
a forth root square dependence as expected for sym- 
metric annular Josephson tunnel junctions. 







